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Abstract. Inelastic effects in electron transport through nano-sized devices are addressed with a method based 
on nonequilibrium Green's functions (NEGF) and perturbation theory to infinite order in the electron-vibration 
coupling. We discuss the numerical implementation which involves an iterative scheme to solve a set of coupled 
non-linear equations for the electronic Green's functions and the self-energies due to vibrations. To illustrate our 
method, we apply it to a one-dimensional single-orbital tight-binding description of the conducting electrons in 
atomic gold wires, and show that this simple model is able to capture most of the essential physics. 
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1. Introduction 

Atomic-size conductors represent the ultimate limit of 
miniaturization, and understanding their properties is 
an important problem in the fields of nanoelectronics 
and molecular electronics. Quantum effects become 
important which leads to a physical behavior funda- 
mentally different from macroscopic devices. One such 
effect is the inelastic scattering of electrons against lat- 
tice vibrations, an issue which is intimately related to 
the important aspects of device heating and stability. 

In this paper we describe a method to calculate the 
inelastic transport properties of such quantum systems 
connected between metallic leads. As a specific exam- 
ple, we here apply it to a simple model for atomic Au 
wires, for which such inelastic signals have recently 
been revealed experimentally (J. 

2. Inelastic transport formalism 

Our starting point is a formal partitioning of the sys- 
tem into a left [L] and a right (i?) lead, and a central 



device region (C), in such a way that the direct cou- 
pling between the leads is negligible. Hence we write 
the electronic Hamiltonian as 

^J^L + yLC + yfc{q) + rR.C + y'PR. (1) 

where M'a is a one-electron description of lead a — 
L/R and YaC the coupling between a and C. The 
central part ^c{<l) is also a one-electron description 
but depends explicitly on a displacement vector q cor- 
responding to mechanical degrees of freedom of the 
underlying atomic structure in this region (within the 
Born-Oppenheimer approximation we assume instan- 
taneous response of the electrons). We are here con- 
cerned with the electronic interaction with (quantized) 
oscillatory motion of the ions. For small vibrational 
amplitudes the g-dependence can be expanded to first 
order along the normal modes A of the structure, i.e. 
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where cj, (Cj,) is the single-electron creation (annihila- 
tion) operator and b\ (b^) the boson creation (anni- 
hilation) operator. The ionic Hamiltonian is just the 
corresponding ensemble of harmonic oscillators 

1. 



(5) 



In the above, the phonon Green's functions Dq^{X, uj) 
are approximated by the noninteracting ones P|- Fi- 
nally, G^'^{lu) are related to Gq^{uj), Sq^(w), and 
Spj^((jj) via the Dyson and Keldysh equations ||2j 



GS(^) + GS(.;)S^(a;)G^(c.), 



(11) 



G^(c.) = G-(u;)[S> + S> + S>](c^)G'^(a;). (12) 



where f^A is the energy quantum associated with A. 

The transport calculation is based on NEGF tech- 
niques 12 ■ For steady state the electrical current 
and the power transfer Pa (per spin) to the device 
from lead a is given by |3] 
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t„(a;) = Tr[S<(a;)G>(c^) 



(7) 

S>(c.)G<(u;)], (8) 



where c/J^ is the electronic number operator of lead 
a. Above we have introduced Green's functions in 
the device region G^(tj) and the lead self-energies 
Sa(ijj) (scattering in/out rates) due to lead a. For 
a shorthand notation these are written as matri- 
ces in the {i^}-basis. For example, the elements in 
G<(w) are the Fourier transforms of G'<(i', t] i/' , t') = 
ih~^{cl,{t')c^{t)). In the hmit of zero couphng 
M^j^, = 0, we can solve exactly for the lead self- 
energies Sq^(cj) and the device Green's functions 
Gq'^{uj) (since this is a single-electron problem). 

Complications arise with a finite coupling, where the 
vibrations mediate an effective electron-electron inter- 
action. To use Eq. and iQ we need the "full" 
Green's functions G'''^((jj). Our approach is the so- 
called self-consistent Born Approximation (SCBA), in 
which the electronic self-energies due to the phonons 
^phi^) are taken to lowest order in the couplings 21 . 
For a system lacking translational invariance |2] 

duj' 
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'D^oi^'^oj- uj')[G<{uj') + G''(t^')]M^ 
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The coupled non-linear equations ©-dl have to 
be solved iteratively subject to some constraints on 
the mode population n\ (appearing in D^(X,lu)). 
We identify two regimes: (i) the externally damped 
limit where the populations are fixed according to the 
Bose distribution ji\ — nsi^x), and (ii) the exter- 
nally undamped limit where the populations vary with 
bias such that no power is dissipated in the device, 
i.e. + Pr = 0. To solve the above we have de- 
veloped an implementation in Python, in which the 
Green's functions and self-energies are sampled on a 
finite energy grid. 

3. Simple model 

As a simple illustration of our method, let us consider 
an infinite one-dimensional single-orbital tight-binding 
chain. We define the central region C to be a piece of it 
with -I- 2 sites to represent the conducting electrons 
in a finite metallic atomic wire. The two semi-infinite 
pieces which surround C can now be considered as left 
and right leads. Ignoring on-site energy and hopping 
beyond nearest neighbors we simply have for C 



^c(g) 



JV+l 

E 



t,,,+i(<?)(c|c,+i+h.c.). (13) 



The hopping amplitudes explicitly depend on the dis- 
placement vector q where the coordinate qi describes 
the displacement of ion i from its equilibrium position. 
As a specific model for the hopping modulation by 
displacement we use the so-called Su-Schrieffer-Heeger 
(SSH) model T in which the hopping parameter is ex- 
panded to first order in the intersite distance 



ti,i+l{q) =t° + t'{qi ~ gi+i). 



(14) 



where i° and t' are site- independent parameters. To 
describe the ions (in a uniform chain where the end 
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sites are fixed in space, qi ~ qN+2 — 0) we include 
only nearest neighbor springs and write 



JV+l 

E 

1=1 



+ \K{qi - qi+if 



(15) 



where m is the ionic mass and K the effective spring 
constant between two neighboring sites. 

Imposing quantization via — iMij, we can 

formulate the linearized electron- vibration interaction 
in terms of the normal mode operators b\ and b 



A' 



-ph 



N N+1 

EE 



i(clq+i-t-h.c.)(6USj, (16) 



and relate the coupling elements to components of the 
normal mode vectors (normalized • = 1) as |3j 



= t'h 



(17) 



It is well established that atomic Au wires have 
one almost perfectly transmitting eigenchannel at the 
Fermi energy (e.g. and references herein). To avoid 
reflection in our model we describe the leads with the 
same electronic parameters as for the wire, leading to 
semi-elliptic band structures of the leads with widths 
4dP . With one electron per site the band is half filled 
and the Fermi energy becomes e^? = 0. Further, we 
take the lead states to be occupied according to Fermi 
distributions np{uj — ^a) where the chemical potentials 
vary as /^l = -1-6^/2 and = —eV/2. With this in- 
formation we essentially have Y^^i^iuj) Pj. The setup 
and the set of normal modes for a particular iV = 6 
atomic wire are shown in Fig. ^ 
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Figure 1: Illustration of the normal modes (longitudinal) of 
a 6-atom wire arranged between two fixed end sites (level- 
broadened due coupling to semi-infinite leads). The open 
circles represent the equilibrium configuration, and the black 
discs a displacement proportional to the normal mode vec- 
tors. The modes are arranged vertically according to the 
mode energy which are also shown to the right of each 
mode vector {K = 2 eV/A ). Note that the highest energy 
mode has alternating bond length (ABL) character. 



4. Numerical results 

Let us now discuss our numerical results for the differ- 
ential conductance calculated with Eq. ijHJl for different 
lengths N and spring constants K. We use the param- 
eter values stated in Tab. ^ which qualitatively yields 
reasonable agreement with the experimental measure- 
ments on atomic Au wires pp. 

Physical quantity Symbol Value 



Bare hopping t'^ 

Hopping modulation t' 

Fermi energy ep 

Atomic mass m 

Spring constant K 

Temperature T 



1.0 eV 
0.6 eV/A 
0.0 eV 
197 a.m.u. 

2.0-8.0 eV/A^ 
4.2 K 



Table 1: Model parameters used for metallic atomic wires. 

The linear energy grid in principle has to cover the 
full bandwidth (FBW) while at the same time it must 
have a resolution fine enough to sample G'^'>{llj) and 
Sa^(w) well. For this model, to resolve the fastest 
variations (caused by the Fermi function) the grid 
point separation should be around 0.4 meV or better 
at a temperature of T = 4.2 K. We find that calcula- 
tions carried out on an interval [— EcutjEcut] converge 
quickly with Scut to those of the FBW. As we show be- 
low for a few representative cases, complete agreement 
is found when Scut = 0.1 eV (which hence are used 
in the calculations presented here). Over this narrow 
range we can further apply the wide band limit (WBL) 
^a{uj) ~ ^ai^ = 0). These simplifications reduce the 
computational load significantly. 
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Figure 2: Differential conductance and its derivative for a 
6-atom wire with different values for the nearest neighbor 
spring constant K in the externally damped limit {nx ~ 0). 
All 6 modes are included in this calculation. 



4 Frederiksen 



The nonlinear conductance versus applied bias 
across a 6-atom wire is shown (i) for the externally 
damped limit in Fig. [3 and (ii) for the externally un- 
damped limit in Fig. |31 It is seen from Fig. |21 that 
the conductance drop essentially happens at one par- 
ticular threshold energy. This energy is found to coin- 
cide with that of the mode with highest vibrational 
energy, i.e. the mode with alternating bond length 
(ABL) character, which can also be designated as the 
dominating one. This mode is further studied in the 
externally undamped limit, Fig. |31 in which a finite 
slope is observed beyond the threshold as well as a lin- 
ear increase in the mode population with bias (heat- 
ing). Generally, both figures show that the conduc- 
tance drop increases while the phonon threshold de- 
creases when the spring constant is lowered. This can 
be interpreted as an effect of straining the wire which 
cause the bonds to weaken. Notice also the agreement 
in both figures between the FBW and the WBL cal- 
culations, shown for the case /-C = 2 eV/A . 



drop also scales with K as (for fixed N), as one 

could speculate from Eq. H17|) . 
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Figure 3: (a) Differential conductance and its derivative for 
a 6-atom wire with different values for the nearest neighbor 
spring constant K in the externally undamped limit. Only 
the dominating mode is included in this calculation, (b) 
Mode occupation n vs. bias voltage. 

With our simple model we can easily handle longer 
wires. In Fig. 0|we show a compilation of the conduc- 
tance drops and the conductance slopes for wires with 
length up to = 40. The individual conductance 
plots all look quantitatively much like those of Fig. [3 
and 121 The important result is that these quantities 
scale linearly with N. If one plots the conductance 
drop against the inverse of mode energy (say, of the 
dominating mode) it is found that the conductance 




(b) 




Figure 4: Compilation of the results obtained for different 
number of atoms in the wire (a) for the externally damped 
limit and (b) for the externally undamped limit. The graphs 
show that the conductance drop and the conductance slope 
beyond threshold scale linearly with the length of the wire. 



5. Conclusions 

In conclusion, we have described a method to calcu- 
late inelastic transport properties of an atomic-sized 
device connected between metallic leads, based on 
NEGF techniques and SCBA for the electron- vibration 
coupling. As a numerical example, we studied a sim- 
ple model for the transport through atomic Au wires. 
With a single-orbital tight-binding description we il- 
lustrated the significance of ABL mode character, and 
were able to explore even very long wires. We further 
discussed the approximations related to a representa- 
tion on a finite energy grid. 

As a final remark, and as we show elsewhere 0, the 
described method is also well suited for a combination 
with full ab initio calculations. The authors thank 
M. Paulsson for many fruitful discussions. 
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